Hadrons as Skyrmions in the presence of isospin chemical potential 
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The stability of the Skyrmion solution in the presence of finite isospin chemical potential /i is 
considered, showing the existence of a critical value /i c = 222.8 MeV where the Skyrmion mass 
vanishes. Using the Hamiltonian formulation, in terms of collective variables, we discuss the behavior 
of different skyrmionic parameters as function of the isospin chemical potential (/x), such as the 
energy density, the isoscalar radius and the isoscalar magnetic radius. We found that the radii start 
to grow very fast for \i > 140 MeV, suggesting the occurrence of a phase transition. 



The skyrmion picture has attracted the attention well as on the temperature. Denning r = eF v r, the mass 
of many authors as a possible way for understanding of the Skyrmion will be given by 
hadronic dynamics and the hadronic phase structure. 
The behavior of hadrons in a media, i.e. taking into 
account temperature and/or density effects, can be ana- 
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lyzed according to this perspective. yr _ r _ ^ j r<\ 

The Skyrme lagrangian is " ~ M=0 ~ Ae^K " 32e 3 F„ [ ' 
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where F n is the pion decay constant and e is a numerical 
parameter. The isospin chemical potential is introduced 
as a covariant derivative of the form |l( 



d I/ U^D v U = d u U-i^[a 3 ,U]g„ . (2) 

The U field matrix, for the static case, can be param- 
eterized in the standard way 



U = Uq = cxp ■ n) = cos£ — i(a ■ h) sm£, (3) 

where a is the sigma matrix vector and h — 1. This 
ansatz has a "Hedgehog" shape obeying the boundary 
conditions 0] 



where Af M= o is the zero chemical potential contribution. 
Notice that the chemical potential terms contribute with 
opposite sign. This implies that the solution will become 
unstable above certain value of fx. In the previous equa- 
tion 
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Assuming a radial profile (£ = £(r), the integrals I2 
and I4 are given by 



£(r) = £(r), h = f, 

e(o)=7r, e(oo) = o. 



(4) 



The mass of the Skyrmion, for static solutions, devel- 
ops a dependence on the Isospin Chemical potential as 
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In order to minimize the mass, we use a variational 
procedure which leads us to the following condition for 
the radial profile 
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FIG. 1: Numerical solution for the profile f(r) for different 
values of /i. (fi = 12.9 (MeV): solid line; fi = 25.8 (MeV): 
dashed line; fi — 38.7 (MeV): dotted line.) 



FIG. 2: The numerical solution for the Skyrmion mass as 
funtion of fi. 



where 
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where /t = ii/(eF^). Equation (JSJ) can be solved numeri- 
cally for the profile £(r), for different values of /i (Figure 
notice that the width of the profile grows with the 
chemical potential. 

In order to obtain the mass of the Skyrmion, the profile 
has to be inserted, numerically, in equation (JSJ. Figure 
[5] shows the chemical potential dependence of the mass. 
The point where the mass vanishes corresponds to the 
critical value fj, c = 222.8 MeV. This value does not de- 
pend on a model for the shape of £(r) . It is a fundamental 
result associated to the Skyrmion picture. 

In , a discussion was presented about the validity of 
the Atiyah-Manton ansatz £| at finite chemical potential. 
It turns out that this procedure is limited only for small 
values of (/i <100 MeV), been in this region quite in 
agreement with the numerical solution. 

The Skyrmion solution we have presented, can be con- 
sidered as the basic or skeleton structure for the hadronic 
states. In order to characterize different hadronic states, 
and following we will introduce collective variables 
A(t), computing then the mass spectrum of the nucleons, 
in an hadronic approach. As usual, the SU (2) collective 
coordinates A(t) are introduced as 



A(t) = a (t)<j + ia(t) ■ a, 
where the a's obey the constraint 

a 2 (t) + a 2 (t) = 1. 
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Introducing and (|lt)(l into Ijl2(l , a rather length com- 
putation leads us to the Lagrangian 



/ Ai 
L = -Ma + 2A di + /i— - 



where A = (27r/3e 3 F 7r ) A, with 
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A = / r sin 2 F 
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In equation 112|l we have defined Ai — gijaj, where Aq — 

03, Ai = —a.2, A2 = a\ and A3 = —ao- is the 

chemical potential dependent mass given in (JSJ). In this 
way we get the Hamiltonian 
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where the canonical momentum is given through a min- 
imal coupling TTi = pi — AXfiAi. The Hamiltonian can be 
expressed as 



U = A(t)U Q A\t). 



(9) 
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The baryonic charge density for the Skyrmion is given by 
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p B = 4irr 2 B°(r) = sin^ F(r)F'(r). (19) 
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Obviously, J °° drps = 1, independently of the shape of 
the skyrmionic profile. The isoscalar mean square radius 
is defined by 
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FIG. 3: The Nucleon static energy as function of the isospin 
chemical potential p. The splitting between neutron (solid 
line) and proton(dotted line) is due to the last term in 11611 . 
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and considering the canonical quantization procedure 

Pi — ► pi = -iS/dai, we get 
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where 7 3 is the third component of the isospin operator 
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Following the usual procedure, we may associate a 
wave function to the Skyrme Hamiltonian. In order to 
identify baryons in this model, these wave functions have 
to be odd, i.e. ip{A) = —ijj(—A). In particular, nucleons 
correspond to linear terms in the a's, whereas the quartet 
of A's are given by cubic terms. 

The energy spectra of nucleons as function of p is 
shown in figure [31 We can see that an energy splitting 
between neutrons and protons is induced. 

The Skyrmion model allows the existence of different 
conserved currents and their respective charges p| . Using 
different charge densities we may define several effective 
radii. The evolution of those radii as function of chemical 
potential provides information about the critical behav- 
ior close to the phase transition, where the Skyrmion is 
no longer stable. Let us first start with the topological 
baryonic current 



This radius seems to be quite stable up to the value of 
p 120 MeV, starting then to grow dramatically. Al- 
though we do not have a formal proof that this radius 
diverges at a certain critical p = p Cl the numerical evi- 
dence supports such claim, as it is shown in figure 

Divergent behavior for several radii, associated to 
different currents, has also been observed in different 
hadronic effective couplings as function of temperature 
in the frame of thermal QCD sum rules [(|. Similar be- 
havior is found for the mean square radius associated to 
the isoscalar magnetic density 
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The divergent behavior of the the effective radii, su 
gests the occurrence of a phase transition, in reference 
the behavior of the nucleons magnetic moments is also 
presented, which also have a divergent behavior. 
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Tr [(U^d„U)(U^d a U)(U^d p U)] . (18) 
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